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Abstract 



A rigorous microscopic calculation of the polarizability of disordered meso- 
scopic particles within the grand canonical ensemble is given in terms of the 
supersymmetry method. The phenomenological result of Gor'kov and Eliash- 
berg is confirmed. Thus the underlying assumptions of their method are jus- 
tified. This encourages application of RMT in the Gor'kov-Eliashberg style 
to more complicated situations. 
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Typeset using REVTgX 



Nearly 30 years ago Gor'kov and Eliashberg (GE) in the pioneering paper |J applied the 
random matrix theory (RMT) and the concept of level correlation to the condensed matter 
physics. They studied the polarizability of a small (mesoscopic) disordered metallic particle. 
The dipole moment of the particle in a weak electric field E(£) is given by the Kubo formula: 

(d i (t))=if (<[*(t) ) di(0]> T )j5i(0* / (!) 

J —CO ^ ' 

Here the square brackets stand for the commutator while the angular ones denote the impu- 
rity (...) and the thermodynamic (• • -}t averages. Introducing a set of eigenstates \k) and 
eigenvalues e& of the Hamiltonian in a particular disorder configuration, one can rewrite Eq. 
(§as 

{dM) = e%(uj) (E ^ : M *i (rMr 3 ) lk ) (2) 
\ M Ek ~ 61 + w + l5 I 

Here (rj)^ is the matrix element of the coordinate: {t^m = (k\fi\l), and h = 1. The 
consequent analysis of Eq. (Q) by GE is based on the following assumptions: 

(i) The product of the matrix elements in Eq.@ can be averaged independently of the 

remaining part of the expression; 

(ii) This averaging can be performed quasiclassically; 

(iii) The remaining summation over the energy levels can be performed with the use of the 
two-level correlation function known from the random matrix theory |2|]. 

This semi-phenomenological calculation inspired a number of attempts to improve it. 
First, Strassler, Rice, and Wyder [|J have pointed out that GE calculated in fact the response 
to the local electric field rather than to the external one. We will return to this point 
later. Then, Devaty and Sievers [|J corrected some minor mistakes contained in the GE 
paper. Shklovskii ||, and later Sivan and Imry |J considered the modification of the GE 
result for the canonical ensemble. Efetov [7[ derived the two-point level correlator for a 
disordered granule by means of the supersymmetry method and showed that it has exactly 
the same form as in RMT, thus having justified the assumption (iii) of GE. Finally, Frahm, 
Miihlschlegel, and Nemeth || attacked the assumption (ii) and tried to improve the GE 
result by introduction of frequency-dependent averaged matrix elements. However, for small 
frequencies uj <C E c with E c = D/L 2 being the Thouless energy (here D and L are the 
diffusion coefficient and the size of the system respectively) this dependence proved to be 
inessential and they have just reproduced the GE result for the matrix elements. 

Below we revisit this problem completely. We restrict ourselves to the case of diffusive 
mesoscopic system (I <C L, with / being the mean free path) and grand canonical ensemble. 
For simplicity we consider only the case of Gaussian Unitary Ensemble (GUE). We present 
a microscopic, well- controlled calculation, based on the supersymmetry method, that does 
not involve GE assumptions (i), (ii), and (iii). We demonstrate that the GE result is correct 
(to some extent, which we also analyze), and hence the assumptions (i) and (ii) are justified. 
Consequently we are able to calculate responses of rather general form, and we discuss the 
conditions of applicability of the GE method. 
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For our purpose the formula (|l|) may be conveniently rewritten in terms of the Matsubara 
Green's functions. In the Fourier representation we have di(u) = aij(w)Ej(u), = aSij, 
and the polarizability a(u) is: 



a[iw n 



— TJ2 J r 1 r 2 dr 1 dr 2 (G(r 1 ,r 2 ,ie m + iu n )G{r 2 ,r 1 ,ie m )) (3) 



The integrand contains a combination rir 2 , which is not translationally invariant. This can 
be cured with the use of the identity 

*i(G ! (ri,r 2 ,e m + a; n )G(r2,r 1 ,e m )) = 0, (4) 

which enables us to replace this combination by a translationally invariant one — (ri — r 2 ) 2 /2. 
After the analytical continuation iu n — > to is performed, Eq. @ is expressed through the 
advanced and retarded Green's functions as follows: 

= — J dedr 1 dr 2 (r a - r 2 ) 2 [n F (e)(G R (r 1 , r 2 , e + u)G R {r 2 , n, e))- 
- n F {e + uj)(G A {r u r 2 , e + cu)G A (r 2 , r ls e)) + 

+ [n F (e + W )-n F (e)](G /? (r 1 ,r 2 ,e + a;)G' A (r 2 ,r 1 e))} (5) 

The products (G R G R ) and (G A G A ) can be calculated in the usual impurity perturbation 
theory (see, e.g. ||). Direct calculation gives: 

J de {n F (e)(G iJ (r 1 , r 2 , e + w)G R (r 2 , n, e)> - n F (e + ou)(G A (r 1 , r 2 , e + u)G A (r 2l n, e))} « 

iv r ( txv \ 2 r 

~ — j expf — -) Ism 2p F r — 2p F r cos 2p F r\ + I exp(— -) cos(2pi?r) (6) 

4r 4 i \Pf^ J I 

Here r = |ri — r 2 |, and v is the density of states. 

At the same time the perturbation theory is not valid for the (G R G A ) term, and we 
involve the supersymmetry technique [f7|,|10| to calculate it. After the standard manipulations 

we get 

(G R G A ){v u v 2 ,e,u) =-J DQ {^(n, r 2 )^(r 2 , n) + «#(n, r0<#(r 2 , r 2 )} e"^ (7) 

Here is the action of the supermatrix sigma model: 

F[Q\ = ^-Jdr Str[D(VQ) 2 + 2i(lu + iS)AQ], (8) 

Q = T~ l hT is a 4x4 supermatrix, A = diag(l, 1, —1, —1), and T belongs to the supercoset 
space C/(l, 1 1 2). The symbol Str denotes the supertrace defined as StrA = AH —^4}/ + ^4 22 — 
Aj 2 f- The upper matrix indices correspond to the retarded-advanced decomposition, while 
the lower indices denote the boson-fermion one. The Green's function g in Eq. ([/]) is the 
solution to the matrix equation: 
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-i(e + ^-Ho)-~(u + W)A + g/2r 



p(ri,r 2 ) = 5(ri - r 2 ) (9) 
We consider now two terms in Eq. (^) separately. The first one yields (cf. ||11|| ): 



(10) 



-(#&&( r i> r a)0» ( r 2, ri)) F « 

„ / [Re G«(n,r a )] 2 - [Im G«(n, r 2 )] 2 (g^ 1 (r 1 )g fc 2 b 2 (r 1 )) F) r « / 
~\0, r>Z 

Here G R is the impurity-averaged retarded Green's functions: 

G R (r u r 2 ) = G A *(r 1; r 2 ) = exp[ip F r - r/2l) 

p F r 

and (■ ■ -)f denotes the averaging with the sigma model action F[Q] (Eq.(PD). 

For relatively low frequencies u <C E c the sigma-model correlators for a closed metallic 
systems can be approximately calculated within the so-called zero-mode approximation J?|. 
This means that only the spatially constant configurations of the field Q(r) are taken into 
account, so that the functional integral over DQ(r) is reduced to an integral over a single 
matrix Q. The latter can be calculated with the use of the technique developed in Refs. 
]7| jrD|jr^] . In particular, the result for the correlator entering rhs of Eq.(|Tl]) reads as 

(QllQlb) = "I " exp — sin — (11) 

with A = (uV)^ 1 being the mean level spacing. 

We will show below that the leading contribution to the polarizability is given by the 
second term in Eq.(^) due to the short-range nature of the first one, Eq. (|10|) . For this reason, 
we do not try to calculate the first term with a better precision and turn to the second one. 
It can be rewritten as 

- (gll{^i)gll{^v 2 )) F = -M 2 (Qll(ri)Qll(r 2 )) F (12) 

In the zero-mode approximation we obtain 

(Qr b (ri)Qtl(r2)} F 

and eventually, 



2iA 



1TL) 



a(u) = "^-Az/ 2 J dr 1 dr 2 (r 1 - r 2/ 
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i.e. the result is frequency independent and real. Thus, in order to obtain the imaginary part 
of a{uj) (which determines the conductivity) and to study the frequency dependence of the 
real part we have to go beyond the zero-mode approximation in treatment of the correlator 
(Qlb(ri)QH(r 2 )) F . For this purpose, we use the method developed in Refs. p^3| , |l4l . The idea 
of this method is to decompose the matrix Q(r) into the constant part Qo (zero mode) and 
the spatially dependent contribution (non-zero modes); then the latter should be integrated 
out. This step is performed perturbatively, with the parameter being g~ x = A/E c <ti 1. The 
next step is to calculate non-perturbatively the integral over the matrix Q . 
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Prior to the calculation it is convenient to rewrite (Q\lQH) in terms of a supertrace: 

(QlK^QlKv,)) = \JDQ Str[ag( ri )6Q(r 2 )] e - F [ Q ] (13) 

Here we have introduced matrices a = diag(l, —1, 0, 0) and b = diag(0, 0, 1, —1). The matrix 
Q can be parameterized as follows: 



Q(r) = T c T 1 A 1 + 



1 + W + -W + -W 6 + - W + ... 



To 
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where the matrix W is block off-diagonal and does not contain a contribution from the 
zero-wave-vector mode: / W(r)dr = 0. Now one has to expand the exponent and the 
pre-exponential factor in powers of W; then the integration over the fast modes can be 
easily performed using the Wick's theorem in the same way as it was suggested in Ref. 
ITS]] . Restricting ourselves to the leading correction to the zero-mode approximation (terms 
quadratic in W), we find: 



(Q%(n)QH(r 2 )) = / dQ Str(aQ &Q 



n(r 1>r2 ) 
4 



Str(aQ )Str(6Q ) -4 



(14) 



71% 

x exp — (a> + i5)Str(AQo) 
The diffusion propagator II is the solution to the diffusion equation 
dn(r ls r 2 ) 



dt 



-DV 2 U(r u r 2 ) = (nu)- 1 5(r 1 -r 2 ) 



(15) 



with appropriate boundary conditions. In particular, if the system is a cube with a size L, 
we obtain 



n( ri ,r 2 ) = £ 



1 



iivV Dq 2 — icu 



q= y{n x , n y ,n z ), m = 0, ±1, ±2, . . . , ^n 2 > 

Since we assume uj E c and the mode q = does not contribute to the sum (|TBp, one can 
neglect u in the denominator. 

Now the integral over Q can be calculated exactly, yielding: 



cosiq^i) cos(g y ?/i) cos^q^) cos(q x x 2 ) cos(q y y 2 ) cos(q z z 2 ); (16) 



(^ 2 (ri)g^(r 2 )) = -2{^ 

1X0J 



I / [ — ^ exp Slu 



n( ri ,r 2 ) 



(17) 



In order to get the final expression for the polarizability we should combine Eqs. (0), 
(|), dp and (|T3). First we stress that the contributions of (G R G R ), (G A G A ), and of the 
term fllOD are short-ranged (localized on distances r = — r 2 | of order of the mean free 
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path I), while that of Eq. ( [LTD is long-ranged. It is easy to check that the contribution of 
the short-range terms to the frequency-dependent part of a(u>) is suppressed for arbitrary 
frequency by the factor il/L) 2 <C 1. Thus, the contribution of the (Q\lQH) term, Eq. flTi 
dominates, and we obtain 



a{uj) 



ie 2 u> 



6n 



\T\V\ 



iA 



-RUJ 



+ 



1+i 



A 

7TU0 



exp — — sm — — 
V A J A 



/(r, - r 2 )^r ld r 2+ 



J (ri - r 2 ) 2 n(r 1 , v 2 )dr 1 dr 2 



(18) 



This expression is valid for any closed or nearly closed |HJ diffusive system irrespective to 
the geometry. For the specified cubic geometry we have 



Re a{u) 



e 2 uL 5 



e 2 (3u t7 



Im a(uS) 



6 7T 6 D 

2e 2 (3u 



Vv 



■Vv 



with 



1 



£ (2* + I) 6 





2tiuj 


(— V 


sin — — 


\7lUjJ 


A 






(— V 


sm — — 


\7TUJ 


A 


7T 6 




~ 960 





(19) 



The result ( |T9"D is in agreement with the GE one (up to the numerical coefficient, caused 
by the difference in geometries). Moreover, the general expression (|18|) is equivalent to the 
corresponding GE expression. This becomes clear when we realize that the time-integrated 
conditional probability introduced by GE (/ dtWr(r' ,t) in their notations) is essentially 
the same as our diffusion propagator n(r,r'). This equivalence confirms the validity of 
the assumptions (i), (ii) and (iii), introduced phenomenologically by GE. It is instructive 
to reformulate these assumptions in terms of the exact single-particle states ipk{?) = \k). 
Then the assumption (i) means that the amplitudes of these states and the corresponding 
eigenvalues may be averaged separately, while the assumption (ii) implies the following 
relation: 



(^(rM(r)^(r'MV)) = V~ 2 U(r, r'), k ? I 



(20) 



As is seen from Eq.(^), the diffusion leads to the long-range correlations of the eigenfunc- 
tions. This effect can not be obtained in the zero- mode approximation (Gaussian ensemble). 

To what extent are the assumptions (i) and (ii) exact? In order to answer this question, 
we have to note that in the same manner as it was done above for the polarizability we are 
able to evaluate any quantity given by 



X(uj) = TJ2 J f(r 1 ,r 2 )dr 1 dr 2 (G(r 1 ,r 2 ,e m + LU n )G(r 2 ,r 1 ,e m )) 



(21) 



with / being arbitrary operator. Alternatively, one can treat Eq. (|21|) by GE method. As 
we have seen, in a particular case / = (ri — r 2 ) 2 the results coincide. It is straightforward 
to see, however, that for / = 1 the exact calculation gives zero due to the identity (01), while 
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the RMT calculation in the GE manner would yield a finite result. The reason is that RMT 
is able to take into account only the long-ranged part of (G R G A ), while other (short-ranged) 
correlators are always lost. The polarizability is determined by the long-distance behavior of 
the integrand, and so the RMT result is correct. However, in cases when the short-distance 
behavior is important, the GE method fails. Summing up, we can say that the GE-style 
calculation gives the correct result only in the case when the dominant contribution to the 
integral Q2ID comes from the large distance range r = |ri — r^l ^> I, which takes place when 
/ grows with r. 

The following comment is appropriate here. Following GE, we have calculated the re- 
sponse to the uniform electric field and neglected effects of screening. The true polarizability 
is however a response to the external field rather than to the local one. Effect of the screening 
on the polarizability was thoroughly studied in Refs. |3],|6|]. Besides, it was noted in Ref. 
that non-uniformity of the local field may be important. We should mention in this context 
that the "local" polarizability Xi defined as the response to the non-uniform electric field 
E(r), 

d(r) = J x(r,r')E(r')dr' 

is given by the expression ( Jl8| ) without the integration over coordinates. This quantity may 
be directly used for the calculation of the true polarizability. 

To conclude, we have provided a rigorous microscopic calculation of the linear response 
expression for the frequency-dependent polarizability of a disordered metallic particle. Our 
result is in agreement with that of Gor'kov and Eliashberg M, thus justifying the under- 
lying assumptions of their derivation. This encourages application of RMT a la Gor'kov- 
Eliashberg to more complicated situations. As a very important example, we mention the 



first steps in non-perturbative treatment of interacting systems p!6yl7| . Since the higher- 
order level correlation functions are involved in that case, an attempt to employ the su- 
persymmetry method meets enormous complications. Therefore, RMT (in spirit of GE 
calculations) is the only available approach to handle these problems. We should stress 
however on the basis of the above analysis that the assumptions analogous to (i), (ii), (iii) 
by GE are not always valid, and in any particular case require a careful investigation. 

The work was supported by the Alexander von Humboldt Foundation (Y. M. B.) and 
SFB195 der Deutschen Forschungsgemeinschaft (A. D. M.). 
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